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The ideal black body fully absorbs all incident rays, that is, all propagating waves created by 
arbitrary sources. The known idealized realization of a black body is the perfectly matched layer 
(PML), widely used in numerical electromagnetics. However, ideal black bodies and PMLs do not 
interact with evanescent fields existing near any finite-size source, and the energy stored in these 
fields cannot be harvested. Here we introduce the concept of the ideal conjugate matched layer 
(CML), which fully absorbs energy of both propagating and evanescent fields of sources acting as 
an ideal sink for electromagnetic energy. Conjugate matched absorbers have exciting application 
potentials, as resonant attractors of electromagnetic energy into the absorber volume. We derive the 
conditions on the constitutive parameters of media which can serve as CML materials, numerically 
study the performance of planar and cylindrical CML and discuss possible realizations of such 
materials as metal-dielectric composites. 


PACS numbers: 81.05.Xj, 78.67.Pt (metamate¬ 
rials), 88.80.ht (wireless transmission), 52.25.Os 
(absorption in plasmas). 

I. INTRODUCTION 

The theoretical concept of the ideal black body which 
completely absorbs all incident ray^ is widely used in 
thermodynamics, optics, and radio engineering. Many 
important applications would benefit from close-to-ideal 
absorbers, including light harvesting, thermal emission, 
stealth, decoupling of antennas and other devices, etc. 
Motivated by these prospects, significant research and 
development efforts have been devoted to realization of 
absorbers which would offer performance close to that 
of the ideal black body, see review paper^^. Most of 
the known results concern planar absorbers, effective for 
a limited range of incidence angles. Omnidirectional re¬ 
alizations using inhomogeneous lossy media have been 
proposed^ and an experimental test at microwave fre¬ 
quencies has been performecP. A theoretical possibil¬ 
ity to perfectly absorb waves of all polarizations com¬ 
ing from all directions is based on the idea of the per¬ 
fectly matched layer (PML) as a layer which produces no 
reflection^^Ki^. All these concepts imply that the body 
fully absorbs all propagating waves incident on its sur¬ 
face, emulating the ideal black body. But are these objects 
truly perfect absorbers or can we perhaps do even better 
and extract energy also from evanescent fields created by 
arbitrary sources? Indeed, the perfectly matched layer 
produces no reflections and fully absorbs power of all 
propagating modes (in the ideal case, realizing the prop¬ 
erties of the perfect black body), but evanescent fields 
of sources are not affected by the presence of the per¬ 
fectly matched layer; thus, no energy is extracted from 
them. Basically, within the known paradigm of uniaxial 


perfectly matched layer^^^, there is a conundrum be¬ 
tween perfect absorption of all propagating waves and 
extraction of power from evanescent fields. If the con¬ 
ditions for perfect absorption of propagating fields are 
satisfied, evanescent modes are not absorbed at all. If 
the parameters of the absorber are modified so that some 
power is extracted from evanescent fields of the source, 
the propagating modes are not any more fully absorbed. 
In other words, evanescent fields of sources contribute 
to absorbed power in lossy media, but within all known 
realization scenarios it appears that if the properties of 
lossy media are chosen so that all the propagating modes 
are perfectly absorbed, evanescent fields do not any more 
interact with the body. 

A typical example is a possibility to extract power from 
near fields of small sources, which has been discussed in 
the context of hyperbolic media electromagnetic^^. Hy¬ 
perbolic media are artificial materials (metamaterials), 
whose permittivity tensor is indefinite, so that the real 
parts of its different diagonal components have the oppo¬ 
site signs. Due to the hyperbolic shape of isofrequency 
contours, evanescent waves exciting a planar interface 
of a hyperbolic material sample are partially converted 
to propagating mode^ within the medium and can be 
effectively absorbed there if the material is loss}^^^^. 
However, lossy hyperbolic bodies are not effective ab¬ 
sorbers for incident propagating waves. Actually, an in¬ 
terface between free space and a hyperbolic medium can 
be matched if the anisotropy axis is titled with respect 
to the interface, but only for one incidence angl^I^. 

In this paper, we study possibilities to create bod¬ 
ies which would act as perfect absorbers for all spatial 
harmonics of incident radiation, both propagating and 
evanescent. Conceptually, such an ideal absorber is able 
to extract infinite power from a finite-size emitter lo¬ 
cated at a finite distance, limited only by the power of 
the source which feeds the emitting object (a generic an- 
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tenna). This study is motivated by the results of theo¬ 
retical studies which prove that a finite-size body, which 
is conjugate impedance matched with respect to its envi¬ 
ronment at all spatial spectrum components, has infinite 
relative absorptivity when compared to the classical ideal 
black bod 3 fflof the same size, e.g. Refs.IHHIl. 

However, to 

the best of our knowledge, only an idealized realization 
of such a conjugate matched absorber as a non-uniform, 
locally isotropic double-negative (DNG) sphere has been 
proposed in RefJ^. In a sense, such object acts as a per¬ 
fect sink for all incident spherical waves characterized by 
arbitrary transverse wave numbers, but its practical real¬ 
ization is challenging because very low levels of material 
losses are required. Here, we show how one can possibly 
create absorbers which are nearly perfectly working with 
both propagating and evanescent incident waves using 
uniaxial DNG materials, both for planar boundaries and 
curved surfaces. In this new scenario, there are prac¬ 
tical possibilities for realizations as well-studied uniax¬ 
ial backward-wave structures (see the discussion in Sec¬ 
tion |v|) . 

Since our proposed structures are not only able to ab¬ 
sorb all propagating waves but can also extract power 
from the evanescent modes, the total extracted power can 
go far beyond what can be attained with know metamate¬ 
rial black holeJ^ and even the ideal black body. In other 
words, we can go beyond “the perfect absorber”. From 
the physical point of view, the conjugate matching con¬ 
dition for evanescent modes is similar to the condition of 
existence of resonant surface modes, which can help to in¬ 
crease exchan ge o f thermal radiation between two closely 
located bodies^®^. However, in the proposed scenario 
it is possible to realize the optimal conjugate matching 
for all modes, both propagating and evanescent, that is, 
create the ultimately effective sink for electromagnetic 
energy from objects located at arbitrary distances. The 
main focus of this paper is on harvesting energy from 
small sources in their near-held zone, but we also discuss 
more general realizations of perfect hnite-size absorbers 
for far-zone sources (plane-wave illumination), that have 
been also proposed in papeil^. 


II. ABSORPTION OF ALL INCIDENT MODES 

A. Eliminating Reflection or Maximizing 
Absorption 


Our starting point is the classical circuit-theory con¬ 
cept of conjugate impedance matching of a load, which is 
the condition for maximizing the power delivered to the 
load from a generator. Referring to Fig. showing an 
arbitrary time-harmonic voltage source V with a com¬ 
plex internal impedance Zq, which is loaded by a com¬ 
plex impedance Zi, we can distinguish two different ap¬ 
proaches to i mpeda nce matching^. The hrst one, illus¬ 
trated in Fig. l(a)[ minimizes (nullihes) reflections from 
the load to the source by equalizing the two impedances: 




(b) 



P 


max 


FIG. 1. Illustration of impedance matching: (a) reflectionless 
matching where the load impedance Zi is equal to the source 
impedance Zq (the power delivered to the load is P), (b) 
conjugate matching where the load impedance Zi is equal to 
the complex conjugate of the source impedance Zq (the power 
delivered to the load is maximized, P = Pmax). 


Zi = Zq. In this case the power absorbed in the load 
equals to (here and in what follows we use rms values for 
the complex amplitudes of the time-harmonic quantities; 
the time dependence is where j = \/^): 


P = 


VP 

4VoP 


Re[Zo]. 


( 1 ) 


The second one, illustrated in 
power 


Fig. Ell maximizes the 


VP 

4Re[Zo] 


that is transferred to the load by satisfying the conju¬ 
gate matching condition: Zi = Zq. Note that in this 
case reflections from the load are in general not zero, be¬ 
cause the matching condition Zi = Zq is satisfied only if 
these impedances are purely real. Obviously, Pmax > P 
if Im[Zo] 7 ^ 0. With the goal to maximize absorption 
in the load, the optimal approach is to ensure that the 
load impedance is conjugate matched with the source 
impedance. 

If the impedances are predominantly reactive (the 
imaginary parts of Zoq are large compared with the 
real parts), the power delivered to the load tends to 
zero when the load resistance tends to zero [eq. Q], 
while for the conjugate matched load it diverges to in¬ 
finity [eq. ^]. In the general electromagnetics context, 
this resonant behavior suggests that conjugate match¬ 
ing wave impedances of two media, in the limit of zero 
losses, unlimited power can be delivered from sources in 
one medium into another medium. This fact was noticed 
in papei^, in studies of maximizing radiative heat flow 
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between two bodies. As a possible approach to realizing 
the conjugate match condition for all modes, in paper^ 
it was proposed to insert dense arrays of thin conducting 
wires into both bodies (orthogonally to the interface), 
which makes the wave impedances of a wide range of 
modes approximately real. This approach is clearly not 
possible if the sources are in free space, and in this paper 
we will study means to ensure all-mode conjugate match¬ 
ing by properly selecting the material parameters of the 
absorbing body. 

Let us consider an infinite planar interface between 
free space and a half-space filled by some lossy mate¬ 
rial, excited by electromagnetic fields created by some 
arbitrary sources in free space, as shown in Fig. |2(a)| 
We can use the concept of conjugate matching to un¬ 
derstand the conditions for maximizing the power ab¬ 
sorbed in the medium. To do that, we expand the in¬ 
cident fields into plane-wave spatial Fourier components 
with the wave vectors k = kou (u = ^Ux + yuy -1- zuz), 
where ko = 27rfy^eoJIo is the free-space wavenumber 
and u is a dimensionless vector. The axis x is di¬ 
rected normally to the interface between free space and 
the absorbing medium and points from free space into 
the medium. To simplify considerations, without loss 
of generality we will consider fields of only one polar¬ 
ization — the TM polarization (magnetic field possesses 
a sole z component and electric field vector lies on the 
xy plane) — and assume that the fields do not depend 
on the coordinate 2 ;. The normal to the interface com¬ 
ponent of the propagation constant in free space reads 

/3o = —jko^u^ — 1 = —jkoSo{uy)^ where Uy G M. When 

\uy\ < 1, the excitation corresponds to propagating waves 
and when \uy\ > 1, we deal with exponentially decaying 
evanescent modes of free space. The corresponding TM 
wave impedance in vacuum (the ratio of the tangential 
to the interface field components) is expressed as 




(b) 


Zo = -jrioM'^y) = -3Vo\Jul - 1, (3) 

where 770 = a/z^o/^o is the free-space wave impedance. 
The branches of all square roots are defined as having 
positive real parts. 

In view of our goal to conjugate match the impedances 
of all modes, we demand the uniaxial symmetry of the 
material tensors of the absorbing medium, with the axis 
directed normally to the interface (along z), to ensure a 
possibility of matching to the isotropic free space. Let us 
write the relative permittivity and permeability tensors 
as 


£rn 

0 

0 


k'rn 

0 

0 

0 

^rt 

0 

, [Mr] — 

0 

k'rt 

0 

0 

0 

^rt 


0 

0 

k'rt 


The propagation constants of transverse magnetic (TM- 
polarized) plane waves in the medium reacP^ Pi = 
—jkos{uy) with 

~ \ ( 3 ) 

V ^rn 


FIG. 2. (a) A uniaxial material sample (region 1, 

£o[£r], yo[yr]) excited by a plane wave (propagating or evanes¬ 
cent) traveling in free space (region 0, eo, Mo) along the direc¬ 
tion of vector u. The corresponding transmission-line model 
representing the vacuum (Pq^Zq) and the sample (/3i,Zi) is 
also depicted, together with the utilized Cartesian coordinate 
system (x,y,z). (b) The variation of the normalized free- 

space TM impedance Zojyo as a function of the normalized 
spectral parameter Uy (the real and imaginary parts). 

where the real part of the square root is positive, ensuring 
decaying waves along the positive direction of the axis z 
(away from the sources). The corresponding impedance 
equals 

Zi = (6) 

^rt 

To understand the conditions for the optimal absorp¬ 
tion of all modes we can use the concept of conjugate 
impedance matching for all plane-wave components of 
the fields. In the configuration of Fig. the load 
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impedance is the wave impedance of a certain mode in the 
medium and the source impedance Zq is the impedance 
of the corresponding (the same Uy) mode in free space. 
The transmission-line model with voltage V (correspond¬ 
ing to the electric field E) and current I (measuring 
the magn etic field H) is also depicted at the bottom of 
Fig.|2( a)| where the two lines with different propagation 
constants (/3o,/3i) and wave impedances (Zo,Zi) repre¬ 
sent the source and the load, respectively. 

The variation of the normalized free-space wave 
impedance (with infinitesimally small losses added to en¬ 
sure the correct choice of the square-root branch) of TM 
waves ^o/^o as a function of the relative tangential com¬ 
ponent of the propagat ion c onstant Uy = (k ■ y)/ko = 
kyjko is shown in Fig. 2(b) (only positive Uy are con¬ 
sidered since all the impedance functions are even with 
respect to Uy). It is apparent that the quantity Zq is 
positive real for \uy\ < 1 (propagating modes) and neg¬ 
ative imaginary for \uy\ > 1 (evanescent modes). With 
the objective to maximize absorption of electromagnetic 
energy by the medium, we will next look for such con¬ 
stituent parameters {Srti that Zi = Z^ for all 

real values of Uy. 


B. Double-Negative Conjugate Matched Layer 

Given the fact that the TM wave impedance of vac¬ 
uum is real for \uy\ < 1, the reflectionless matching 
is identical to conjugate matching for the propagating 
waves. Therefore, for that part of the incidence spec¬ 
trum, we are searching for a uniaxial medium that does 
not reflect incoming traveling fields. Such a medium is 
the famous uniaxial perfectly matched layer (PML) witlP 

^rt ~ k'rt ~ • (^) 

£rn 

Indeed, substitution of 0 into 0 and comparison with 
® shows that Zi = Zq for any value of Uy. However, 
the fact that a half space with these characteristics ex¬ 
hibits no reflections also for evanescent modes, namely 
that Zi = Zq also for \uy\ > 1, does not fit well with the 
maximal absorption goal. 

It is known that a planar interface with an isotropic 
Veselago medium {cr = yUr = ~1) is conjugate matched 
with free space for all modeJ^, which is the enabling 
property for the perfect lens operatioiP^. However, this 
material is lossless for both propagating and evanescent 
waves. If we introduce small losses in order to extract 
energy from evanescent waves, the perfect matching for 
propagating modes is destroyed. The power generated by 
evanescent fields can be increased by making losses small 
[eq. ©]. but at the same time absorption of propagating 
modes gets smaller and smaller (the imaginary part of the 
propagation constant tends to zero). To overcome this 
deficiency, we propose to use double-negative materials 
whose parameters satisfy 0> but Re[6^t] = Re[yU^t] < 0- 
Obviously, the real part of the normal component of the 




(b) 

FIG. 3. The variation of: (a) the real parts and (b) the 
imaginary parts of the normalized TM impedance Zi/rjo for 
a uniaxial medium as a function of the normalized spectral 
parameter Uy. We examine three cases (routes) based on the 
PML concept: (i) normally matched double-negative (DNG) 
isotropic medium with Em = Srt = Mm (isotropic route), (ii) 
DNG lossy PML-type medium with a lossless normal compo¬ 
nent, where: Em = Re[l/£rt] = Re[l/Mrt], so that none of 
the three constituent parameters (ert, is active (non¬ 

active route) and (iii) DNG lossy PML-type medium with 
Ern = 1/^rt = (activc Ern, PML routc). In all the cases 

Ert — l^rn — 1 4 0 - 3 . 


permittivity is also negative: Re[5r^] = Ile[l/ert] < 0. 
Using such materials, conjugate matching for all propa¬ 
gating modes can be realized also in lossy configurations, 
where all propagating plane waves quickly decay inside 
the absorbing medium. This property allows us to expect 
that we can effectively extract power from both propagat¬ 
ing and evanescent waves. We will call such structures 
(double-negative) conjugate matched layers (GML). 

It is needless to say that if the constituent parame¬ 
ters are purely real, no absorption takes place. Accord¬ 
ingly, we should add some losses to the transversal com¬ 
ponents; in this situation, all propagating waves expo¬ 
nentially decay and give up their energy to the medium. 
However, in this case the normal permittivity compo¬ 
nents Ern = 1/^rt = ^/hrt ^rc inevitably active. Indeed, 
the sign of the imaginary parts of Em and (im is the 
opposite to that of the imaginary parts of Ert and {irti 
as the material parameters comply with the PML con¬ 
ditions 0 . Thus, the material parameters of the ideal 
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double-negative conjugate matched layer do not have any 
active components only in limiting cases when losses tend 
to zero or the real parts of the transverse components are 
very large. 

Let us next show that double-negative materials whose 
parameters satisfy 0 are indeed conjugate matched to 
free space for all modes. For propagating waves {\uy \ < 1) 

we have s{uy) = yj— 1) = which 

corresponds to the normal component of the propaga¬ 
tion constant f3i = —jkos{uy) having a negative real 
part (propagating backward waves, as in any double¬ 
negative material). The chosen branch Re(^) > 0 corre¬ 
sponds to the proper direction of the power flow from the 
sources into the absorbing medium. The corresponding 

impedance reads = —jr]os{uy)/ert = — 1 = 

Zq = Zq, because this is a real number. For evanes¬ 
cent modes {\uy\ > 1) we get s{uy) = = 

—Srt'^Uy — 1, where the sign is chosen so that the fields 

decay away from the sources in free space. The cor¬ 
responding impedance becomes = —jr]os{uy)/ert = 


JVo 


ul-l = Zo*. 


In Figs, [^we show the v ariation of the real (Fig. 3(a)| ) 
and imaginary (Fi g^3(b) ) parts of the normalized wave 
impedance Zijrjo Q as functions of Uy for materials 
with the parameters approximating the ideal values of 
the CML layer. We examine three possible routes to¬ 
wards CML: (i) An isotropic DNG material with Sm = 
Srt = iirn^ which has the CML properties in the limit 
Crn = ^rt = l^rn ^ “1; fhis approach is labeled 
as the “isotropic route”, (ii) A lossy uniaxial double¬ 
negative material with the lossless normal component of 
the permittivity, satisfying Sm = Re [l/srt] = Re [1/ 
which formally has the ideal CML properties in the limit 
Re[5^t] = Re[/irt] ^ —eo; this approach is labeled as 
the “nonactive route”. Finally, (hi) a double-negative 
lossy CML medium with Sm — ^/^rt = fhis ap¬ 

proach is labeled as the “PML route”. The impedances 
are compared with the ideal complex conjugate wave 
impedance Zq of free space: The equality of the two 
quantities guarantees maximal absorption. 

It is clear that the choice of negative real parts for 
{Srt^ l^rti^rn) Is succcssful in terms of ensuring the de¬ 
sired sign of Im[Zi]; one should take into account that 
for the usual PML with positive real constituent param¬ 
eters (double-positive, DPS), we obtain = Zq for 


all real Uy and thus Im[Zi] < 0, as seen in Fig. 2(b) 


It is also apparent that for both real and imaginary 
parts, the uniaxial material with lossless £rn performs 
better than the isotropic material. Most importantly, 
the double-negative lossy CML corresponds to exact con¬ 
jugate matching^ the fact that verifies our initial as¬ 
sumption. Note that such a selection of a “PML” with 
Re[5rt] < 0 and lm[€rt] < 0 changes (compared to the 
normal DPS PML) only the wave impedance of evanes¬ 
cent modes since Zi is purely real along the propagating 


spectrum. 

One can say that we misapply the concept of PML 
since the choice of double-negative components inflicts 
strong reflections {Zi = Zq) for the evanescent modes, 
namely for \uy\ > 1, contrary to the usual reflectionless 
matching property of DPS PML {Zi = Zq). However, 
within our paradigm, all propagating waves are fully ab¬ 
sorbed, and evanescent fields, in the limit of ideal con¬ 
jugate matching, generate infinitely strong fields in the 
medium and deliver infinite power to the absorber in the 
limit of Re[Zi] ^ 0 [Eq. 0], in the assumption that the 
illuminating antennas are fed by ideal voltage or current 
sources capable to supply unlimited power. 



FIG. 4. A grounded uniaxial slab of thickness L is excited by 
an electric dipole line source forming an angle 0 with the hor¬ 
izontal axis, placed at a close distance g from the boundary. 
The corresponding transmission-line model is also depicted. 


III. INFINITE PLANAR CONFIGURATION 
A. Incident Field and Absorbed Power 

To test and demonstrate the absorbing efficiency of 
double-negative conjugate matched layers, we consider 
the planar structure depicted in Fig.[^ A grounded slab 
of thickness L (region 1) filled with a uniaxial medium 
{srt^ l^rt^ ^rn) IS cxcitcd by an electric dipole line sourc^^ 
located at the point (x, y) = (—^, 0) in the vacuum region 
(region 0). The dipole moment per unit length of the line 
is denoted as p/. The dipoles form an angle 0 with the 
horizontal axis x. The distance between the slab and the 
source is much smaller than the free-space wavelength 
(in order for the evanescent waves not to vanish before 
reaching the sample), and is denoted by To study 
absorption of evanescent fields, we choose such a source 
since its evanescent spectral components are much more 
pronounced compared to that of more typical cylindrical 
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wave sources such as lines of electric or magnetic cur¬ 
rents!^ The problem is effectively two-dimensional (2D), 
due to ^-independence, and the magnetic field has a sin¬ 
gle z component. 

We write the solutions in spectral domain, as a plane- 
wave deco mposition in the form exp(—The expres¬ 
sion ( |A.4[ ) of the Fourier-transformed incident magnetic 
field ?io,inc (^5 ky) = z?^o,inc (^5 ^y) obtained in Appendix 
can be evaluated at x = 0 as follows: 


^0,inc(^ — 0? ky^ 


2 


~ k 

cos 0 — sin 0 
_Po 


( 8 ) 


where pi = |p^|. This field excites the uniaxial slab lo¬ 
cated at X > 0. We note that identical plane-wave field 
components are generated by sheets of surface electric 
current = — 2 yHo,inc (05 placed right in front of 
the grounded slab. Therefore, the excitation of the slab 


by the line of electric dipoles (Fig.[^ top) can be equiva¬ 
lently represented by the transmission-line model shown 
at the bottom of Fig. with the external current source 
Js inserted at x = 0. Using this circuit model, the power 
delivered to the slab by each incident spatial harmonic is 
straightforwardly calculated as: 


P{ky) = 


m0,inc{0,ky)\^\Z0{ky)\^ 
\Z0{ky) + Zi^{ky)\^ 


Re[Zin{ky)]. 


(9) 


As above, by Zq we denote the free-space wave 
impedance, and Zi^ is the input impedance of the metal- 
backed uniaxial slab evaluated from the corresponding 
transmission-line model as Zi^ = jf’Zi tan(/dil/). Finally, 
the total power absorbed in the slab is found by integra¬ 
tion over all the wavenumber spectrum: 


1 

P = ^ J P{ky)dky 


, ,2^2 r-\-oo 

I (^y)9 

27r y_oo 


ky 

Mky) 


COS 0 — sin 0 


2 


Re[Zir,{ky)] 

-I Zinjky) 
^ Zoiky) 


dh 




( 10 ) 


Here it is denoted = — Im(/do) (for propagating waves, 
;do = 0, and for evanescent waves, = \ k‘^ — k^). 


B. Important Limiting Cases 


Let us now consider some limiting cases under the sim¬ 
plifying assumption of a half space emulating an electri¬ 
cally thick absorbing slab (valid when I3'{L ^ 1, where 
Pi = —Im[/3i]). In this case, the input impedance 
of the metal-backed uniaxial slab approaches its wave 
impedance: Zi^ = Zi. In our study, the material param¬ 
eters of the uniaxial layer are chosen so that the conju¬ 
gate matching condition is closely approached: Zi = Zq. 
When \ky\ < ko (propagating waves), Zq = r]oPo/ko is 
purely real and, thus, = Zq. Therefore, the power 
delivered by the propagating modes to the conjugate- 
matched layer can be calculated as: 



rko 

k 

cos 6* — sin 6* 


f-ko 

/3o 


^dky 
ko ^ 


Pquj^ 

16 


( 11 ) 

Note that this result is independent of the angle 6 be¬ 
cause (11) is nothing more than 1/2 of the total power 


emitted by a line of dipoles in unbounded free space. This 
power also equals the total amount of power absorbed by 
an ideal PML layer (or, equivalently, by a perfect black 
body) in the structure shown in Fig. 

When \ky\ > ko (evanescent waves), Zq is purely 
imaginary and, thus, fulfilling the conjugate matching 
condition results also in almost purely imaginary Zi: 
Z\ = —jlni[Zo]. In this case, the denominator 


of (10) approaches zero. Therefore, in order to analyze 


this limiting case we must consider a particular model 
for the material parameters of the uniaxial layer. Here 
we select the PML route with Srt = k^rt = a — jb and 
Ern = ^l^rt — where real 5 is such that \6\ 1. 

Performing the Taylor expansion with respect to small 
6 under the aforementioned assumption \ky\ > ko and 
keeping only the first-order term, we obtain the follow¬ 
ing approximation for Zi: 


Zl{ky) 


VoPi 

Srfko 


^ sgn{a)Zo{ky) 


l^jS 


{a - jh)kl 

2{kl-kl) 


( 12 ) 

Separating the real and i maginary parts and taking into 
account that ^/(a — jb^ = sgn(a)(a — j6), we come to 
the following formula: 


Ile[Zi{ky)] 


\a\k. 


\l^Zl{ky)/Z0{ky)\^ 


2koJkl - kl 


a6k^, 


1 2 


‘^{^y ^o) 


1 + sgn(a) + 6sgn(a) 


(5/c2 


2(^y ^o) 


0 . 


(13) 


Now, we can substitute (13) into (10) and calculate the integral over the evanescent part of the spectrum {\ky \ > ko). 
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This calculation results in the following approximate expression for the power absorbed from the evanescent modes; 


Hp r 

°Jko 


8\a 

kn 


[k: 


kl sin^ 


(h‘^ _ ^2 


3/2 


/le-29v^ 


2_U2 

n, 


[1 + sgn(a)]^ + (52 


kl\e^t\ 


2 ) 


0 . 


2(fe2-fcg) 


(14) 


In the denominator of the integrand of (14) we have 


dropped the term proportional to 6, because in the DNG 
case (a < 0) it is multiplied by 1 + sgn(a) = 0 and in this 
case the first nonvanishing term is the second-order term 
in 5 [this term is present in ([T^], while in the DPS case 
(a > 0) the main term is [1 + sgn(a)]^ = 4 ^ 5 ^ 5^. 
Equation ( p!4| ) demonstrates the dramatic influence of the 
sign of the real part of Srt = l^rt- In particular, for the 
CML case when a is negative, the denominator of (14) 
is approximately proportional to and the following 
formula results from (1T4|: 


Po 


32|a| 

TS\Srt\‘ 


r-\-oo 

L ' 


(1 




dso, 
(15) 

where a change of the integration variable sq = 
yjky — k^/ko = — 1 has been made. We notice 

that the absorbed power is proportional to 1/(5, which 
means that for small S we have huge magnitudes either 
of absorbed {S > 0) or emitted {6 < 0) powers. Such a 
limiting regime is explained by the singularity of for 
Re[^i] —y 0. 

Expression (14) can be applied also for the nonac¬ 
tive route when Sm = Ee(l/5pt)- Indeed, when 6 = 
6/(a^ + 6^) ^ 1, the imaginary part of Sm vanishes 
and the normal component of the permittivity dyadic 
acquires the value Sm = + b‘^) = Re(l/5pt)- Re¬ 

spectively, Eq. (p^ becomes 
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In this case, the absorbed power increases with the de¬ 
crease of b and the increase of \a\. 

On the other hand, for the conventional DPS PML 
case, expression (14) takes the following form: 
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(17) 

In this case, when (5^0 (approaching the ideal 
impedance matching condition), the relative power de¬ 
livered to the slab by the evanescent waves tends to zero. 
This result confirms that within the know PML scenario, 
the goal of full absorption of propagating waves is incom¬ 
patible with the goal of extraction power from evanescent 
fields. 

Note that the above approximate formula makes sense 
only when 0 = 7r/2, because only in this case the inte¬ 
gral 0 converges. The divergence of the integral ([T^ 


I 

at any other values of 0 is an artifact of the Taylor ex¬ 
pansion approximation (12). This approximation is not 
applicable when ky approaches ko (i.e., when sq ^ 0). 
However, direct evaluation of the integrand in (p^ shows 


no singularity at ky = ko, even when 9 ^ 7r/2. There¬ 
fore, when 0 ^ 7r/2, the lower limit of the integration in 
Eq. (17) must be replaced by a ^ because at 


this point the first term of the Taylor expansion (12) is 

I \ 

about unity: 2 {kJ-^ 2 ^^ b = 1. In this case, Eq. (17) quali¬ 
tatively agrees with the accurate result ( p!Q| ) (wEen only 
the contribution of the domain \ky\ > ko is considered), 
with the dominant contribution to the integral coming 
from a narrow region in the vicinity of the lower limit. 
Therefore, when 0 ^ n j2 and kog ^ 1, 
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l^rt I 


i.e., the power delivered by the evanescent modes in this 
case is roughly proportional to VS and also vanishes when 
(5^0. 


C. Performance of the Proposed Structure 

Here we characterize the proposed conjugate matched 
layers qualitatively, studying the ratio of the total ab¬ 
sorbed power and the power which can be harvested only 
from the propagating part of the spatial spectrum of the 
excitation in the layer of the infinite thickness (11). Basi¬ 


cally, this normalized parameter is the ratio of the power 
absorbed by our structures and the power absorbed in the 
classical ideal black body at the same position and excited 
by the same source. In particular, the output quantity 
is the logarithm (base 10) of this ratio log(P/Po) since 
we are expecting huge magnitude variations of |P|, as 
discussed in the previous section. 

In Eigs. we represent the absorption enhancement 
parameter log(P/Po) on the complex permittivity map 
(a, 6) with Crt = Prt = a — jb for the three considered 
routes. In the isotropic scenario (Eig. |5(a')] ), the overall 
absorption is small while a local maximum is exhibited 
along the line a = — 1, which corresponds to the Vese- 
lago materiaP^. This result is expected, because in the 
absence of losses this mediu m is co njugate matched for 
all modes at this single poinfP^^. As far as the non¬ 
active route (Eig. 5(b)[ ) is concerned, we observe much 
better performance compared to the isotropic route: the 
power extracted from the source is large in a wide range 











































FIG. 5. Contour plots of the quantity log(P/Po) on the 
complex plane of the transversal constituent relative pa¬ 
rameters (a = Re[£rt] = Re[firt],b = —lm[£rt] — 
for the three considered routes. (a) The isotropic route 
{srn = Srt = l^irt = a — jb). (b) The nonactive route 
{ern = Re [l/ert] = Re [l/^J-rt] = a/{a^ + (c) The PML 

route {ern = llsrt - j5 = 1/Urt - jS = l/{a - jh) - j5). 

Plot parameters: /cqT = 20, k^g = 0.5, 0 = 90°, 5 = 0.001. 


of material parameters, and it is monotonically growing 
when \a\ increases along the n egative semi-axis (a < 0). 
Most importantly, in Fig. |5(c){ where the PML case is ex¬ 
amined, we clearly notice a huge change in the absorbed 
power P from double-positive (a > 0) to double-negative 
(a < 0) cases, as expected from ( p^ . In fact, extremely 
high enhancement of absorption is achieved for the CML 
case (P = lO^Po) which means that the structure is suck¬ 
ing all the power (for every single mode, propagated or 
evanescent) from the source. 

In Figs. and ^ we show the variation of log (P/Pq) 
with respect to the electrical thickness of the grounded 
slab koL and the loss factors b of its relative constituent 
parameters for the three routes (isotropic, nonactive, and 
lossy PML) for different signs of a (Figs. [^for a > 0 and 
Figs. for a < 0). In th e isotropic scenario of the double¬ 
positive case (Fig. |6(a)'] ), we observe resonances when k^L 
is moderate and the reflected fields from x = L are quite 
strong to interfere with the incident waves, while the ab¬ 
sorption P is positively related to the loss param eter b . 
When it comes to the DPS nonactive scenario (Fig. |6(b)"] ), 
we have very small absorption for 6 ^ 0 (as physically 
anticipated), but there is an optimal range of 5, for each 
fixed koL across which the absorbed power gets maxi¬ 
mized. It is also remarkable that the DPS PML config¬ 
uration (Fig. 6(c)) performs worse than the other two 
in terms of absorption; furthermore, note that P is not 
substantially dependent on losses 6, as is clear from 


FIG. 6. Gontour plots of the quantity log(P/Po) with 
respect to the electrical thickness of the grounded 
slab koL and the losses of the transversal constituent 
parameters b = —lin[£rt] = —lin[firt] for the three 
considered DPS routes. (a) The isotropic route 
{srn = £rt = l^rt = a — jb). (b) The nonactive route 
[sm = Re [l/srt] = Re [l/grt] = Re [1/(a - jb)] = a/(a^ + 5^)) 
(c) The PML route (sm = 1/^rt — jb = 

1 / grt — jS = 1/(a — jb) — jS). Plot parameters: 

a = Re[£rt] = Re[Mrt] = 2, kog = 0.5, 0 = 90°, S = 0.001. 


(18). As far as the double-negative configurations are 
concerned, we remark that the behavi or of the struc¬ 
ture following the isotropic route (Fig. |7(a)| ) is similar 
to that o f the corresponding double-positive configura¬ 
tion (Fig. |6(a')]) . On the contrary, the nonactive realiza¬ 
tion (Fig. |7(bJ| ) does well for small losses b and its ab¬ 
sorbing efficiency deterior ates with increasing b. Finally, 
the CML structure (Fig. |7(c)| ) absorbs, on the average, 
extremely high power P, as predicted by (15), and the 
harvested power is practically not fluctuating with the 
electrical thickness k^L^ because the evanescent mode ab¬ 
sorption is due to resonant surface modes. 

In Figs.j^we represent the logarithm of the relative ab¬ 
sorbed power log (P/Po) with respect to the electrical dis¬ 
tance from the source to the boundary k^g. All the three 
scenarios are considered for do uble-positive (Fig. 8(a)) 
and double-negative (Fig. |8(b)[ ) configurations. One can 
notice that the DPS PML absorbs power Pq only from 
propagating waves since a > 0 and (5 ^ 0, as sug¬ 
gested by (17), (18); in other words, PML is outper¬ 
formed by the other two scenarios in the double-positive 
paradigm. Again, we remark a similar response of the 
isotropic structure, regardless of the source location k^g 
(the blue curves in the two figures), while the nonactive 
configuration performs always in-between the other two 
scenarios (which are switching positions when changing 
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FIG. 7. Contour plots of the quantity log(P/Po) with 
respect to the electrical thickness of the grounded 
slab k{)L and the losses of the transversal constituent 
parameters b = —lm[srt] = —lin[firt] for the three 
considered DNG routes. (a) The isotropic route 
{srn = Srt = /J^rt = a — jb). (b) Thc uonactivc route 
[sm = Re [l/srt] = Re [l/^rt] = Re [l/(a - jb)] = a/(a^ + b‘^)) 
(c) The PML route = l/srt — jb = 

— jb = l/(a — jb) — jb). Plot parameters: 

a = Re[e:rt] = Re[/Xrt] = —2, kog = 0.5, 0 = 90°, 

(5 = 0.001. 




the sign of a). In Fig. |8(b)| we observe huge absorp¬ 
tion for the DNG CML case at moderate distances to 
the source (P = IO^Pq), as expected from (15). Finally, 
when kog Too, according to all the scenarios, the de¬ 
vice tends to absorb power P = Pq which indicates that 
only propagating waves survive and excite the boundary. 
The magnitude is the same since, due to the position of 
the source {0 = 90°), the illumination is almost normally 
incident on the grounded slab and only the transver¬ 
sal relative parameters (the same for all the scenarios, 
Srt = jJirt) 9.re activated. 

To conclude this section we note that although we have 
considered only one (TM) polarization, all the results are 
general and hold also for the TE polarization, which can 
be studied in the same way. 


D. Comparison with Hyperbolic Media 

Probably the only known alternative approach to en¬ 
hance absorption of evanescent fields in large lossy bodies 
is the use of hyperbolic medisP^K^. Due to the trans¬ 
formation of evanescent waves to propagating ones in 
hyperbolic media, subs tantial absorption enhancement 
has been observe d Therefore, it would be mean- 


FIG. 8. The quantity log (P/Po) as function of the electrical 
distance of the source kog from the boundary for the three 
considered routes, (a) Double-positive cases a = Re[£rt] = 
Re[/irt] = 2. (b) Double-negative cases a = Re[£rt] = 
Re[//rt] = —2. Plot parameters: b = —lin[ert] = —Im[/irt] = 
0.5, koL = 20, 6> = 90°, b = 0.001. 


ingful to compare our results with those obtained when 
utilizing hyperbolic materials. In hyperbolic media, the 
real parts of different eigenvalues of the permittivity ten¬ 
sor have opposite signs. Thus, we consider the same 
structure of Fig. with Srt = cl — iirt = 1 and 
^rn = for a > 0 and small 6 > 0. We assume 


that Re[£^t]Re[^rn] = “1? because the maximal absorp¬ 
tion is achieved under this assumptiorP^. In Figs. we 
present the relative wave impedance Zijrjo (the real and 
imaginary parts) of the hyperbolic material as functions 
of the normalized spectral parameter Uy for several a. 
It is apparent that along the entire Uy axis, the differ¬ 
ences from the ideal ^o/^o are substantial, which means 
that the absorption would be dramatically smaller than 
in our conjugate-matching scenarios. This expectation is 
confirmed by the results shown in Figs.[^ where the log¬ 
arithm of the absorbed power log (P/Pq) is presented as a 
function of the real part of the permittivity a for va rious 
loss factors b. In the hyperbolic scenario (Fig. |lQ(a) ), the 
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(a) 



(b) 

FIG. 9. The variation of: (a) the real parts and (b) the 
imaginary parts of the normalized TM impedance Z\ /770 for a 
uniaxial hyperbolic medium with Srt — ol — jh, firt = 1, Sm = 
— Ija — jb as function of the normalized spectral parameter 
Uy for various a. The ideal for conjugate matching impedance 
Zq /r]o is also shown. 


power absorption is, on the average, two orders of mag¬ 
nitude smaller as compar ed to t he absorbed powe r along 
the DNG nonactive (Fig. |lQ(b)] ) or the PML (Fig. |lQ(^ ) 
routes. Furthe rmore, one can p oint out that for large \a\ 
(a < 0 for Figs. |lQ(1^ and [TQ(^ the nonactive route gives 
better outcomes than the PML route. This property fol¬ 
lows from the fact that we have assumed a constant ad¬ 
ditional loss factor S for the PML route. In particular, 
for \a\ +OC, the normal permittivity component of 
the nonactive material tends to zero (which is the ideal 
value in the PML route case, corresponding to conjugate 
matching), but for the PML route, we have Sm ~j^- 
In other words, when |a| is increasing without limit while 
6 and b remain fixed, the CML response to evanescent 
excitations appears more lossy than the response of our 
nonactive device. 


Hyperbolic Route DNG Nonaetive Route 




DNG PML Route (CML) 



FIG. 10. The quantity log (P/Pq) as function of the transver¬ 
sal permittivity a for several losses b, evaluated for: (a) the 
hyperbolic case Srt = a — jb^firt = l,£rn = —Ija — jb with 
a > 0. (b) the nonactive route case Srt — IJ^rt = a — jb, Sm — 
Re [l/(a — jb)] = a/(a^ + b‘^) with a < 0. (c) the PML route 
case Srt = Mrt = a — jb,£rn = l/{a — jb) — jS with a < 0. Plot 
parameters: koL = 10, kog^O.b, 0 = 90°, 6 = 0.001. 


IV. FINITE CYLINDRICAL CONFIGURATION 

A. Conjugate Matched Cylinder via Coordinate 
Transfer mat ion 

Let us next study possibilities for the use of the double¬ 
negative conjugate matched materials to create ideal ab¬ 
sorbing bodies of finite sizes (in cross sections normal 
to z). In contrast to an infinite planar surface, stud¬ 
ied above, finite-sized bodies can in principle absorb in¬ 
finite power carried by a single propa gating plane wave, 
generated by sources at infinityl^®^. A conceptual ex¬ 
ample of a spherical body having these properties has 
been describee?^, where the sphere was filled with a lo¬ 
cally isotropic low-loss double-negative material. As dis¬ 
cussed above, double-negative materials can be ideally 
conjugate matched to free space, but only in the limit of 
negligible loss factors. Here we study a cylindrical config¬ 
uration and explore additional design possibilities offered 
by double-negative perfectly matched layers, which have 
more free parameters to tune. 

With the goal of finding the material parameters 
of a cylindrical CML, we adopt an already proposed 
approaefP^ and consider first a complementary infinitely 
long cyli ndrica l shell filling the space from p = r to p = oc 
(see Fig. msi where the used cylindrical coordinate sys¬ 
tem (p, (p, z) is also defined). We first assume that the 
shell is filled with a medium characterized by the follow- 
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Mo) 






(a) (b) 


FIG. 11. (a) Complementary cylindrical infinitely thick man¬ 
tle with the internal boundary p = r hlled with the con¬ 
ventional PML relative constituent parameters (£o[£r], Mo[Mr]) 
located in the unprimed coordinate system (b) 

Cylindrical core with the external boundary p' — r hlled 
with a material having transformed constituent parameters 
(5o[£r]? Mo[Mr]) of (20), locatcd in the primed coordinate sys¬ 
tem {p,ip\z'). 


ing, more general (compared to the PML case) permittiv¬ 
ity and permeability dyadics in cylindrical coordinates: 

[^r] ~ [/^r] ~ '^r2:) • 

In order to obtain the corresponding material param¬ 
eters for a cylinder with a hnite radius r, we map the 
outer space of the cylindrical volume (p > r) into the in¬ 
ternal v acuum hole (p < r) and vice versa, as shown 
in Fig. |ll(b)[ This mapping can be done using the 
concept of transformation optic^^^^, and it has been 
already performed in the analogous case of isotropic me¬ 
dia and the spherical shape of the conjugate matched 
bod}®!. More specihcally, we can use the following coor¬ 
dinate transformation from the unprimed (p, p, z) to the 
primed {p'^ip'^z') coordinates: 


^n+l 

p' = V>' = ‘P, z' = -z, (19) 

P^ 


where n can be any positive real number. Note that the 
sign of the ^(-component changes in the transformation. 
As p' is inversely proportional to p, the corresponding 
coordinate axes have opposite directions. Therefore, in 
order to map a right-handed coordinate system to a right- 
handed coordinate system, we also change the direction 
of either the p or z axes. 

The material parameters transform according to the 
well-known coordinate transformation equations^ 


^ A • [gr] • ^ A ■ [llr] ■ A^ 

det(A) ’ det(A) 


( 20 ) 


where denotes the transpose of a matrix, [sr] and [p^] 
are the original material parameter matrices (i.e., the 
material parameters of the infinite mantle at p > r), [e'^] 
and [p^] are the transformed material parameter matrices 
(i.e., the material parameters of the finite core, p < r), 
and A is the Jacobian matrix of the transformation. The 


Jacobian matrix in this particular case of (19) is given by: 


A = 


1 

1 dp' 
p dp 
p_ d^ 

_ 1 


-n 1 

(5)”‘ 

0 

> 

0 

F dp 
dz' 

_ dp 

p dp 

1 dz' 

P 

F dz 
dz' 

dz J 



0 

0 

p_ 

p 

0 

0 

-1 _ 


( 21 ) 

By inserting the above expressions into (20), the trans¬ 


formed relative material parameters for the core can be 
written in cylindrical coordinates as: 


[4] = [4] = 


ns 


0 

Q Srcp 


rp 


0 


0 ^ A 




\ 2 -|- 2 /TL 


( 22 ) 


Furthermore, if we assume that the initial cylindrical 
shell has PML characteristics, that is, its relative con¬ 
stituent parameters are defined as Sr(p = Srz = 1/^rp = 
a — jh (still assuming [sj] = [pr]), we can write the trans¬ 
formed parameters as 


[4] = [P'r] = 


a-jb 

0 

0 


0 

a-jb 

n 

0 


a-jb 


{7) 


2+2/n 


(23) 


It should be noted that we can derive the same material 
parameters if we, alternatively, demand cp' = —cp instead 
of z' = —z in the initial coordinate transformation. At 
this stage we are free to choose any real values for a, 6 G 
M and any positive value for n; however, in this study 
we limit the analysis to the case n = 1. Furthermore, 
we confine ourselves to the TM case with no harm to 
generality. 

As we did in the case of an infinite planar slab, we 
consider three different routes approaching the ideal con¬ 
jugate matched cylinder, namely the isotropic, nonactive 
and PML routes. Here, the isotropic route corresponds 
to isotropic shell parameters, that is, Sr = Pr = a — jh. 
Clearly, according to (22) this corresponds to trans¬ 


formed parameters that are not isotropic. However, we 
will use the term “isotropic” also in this case to main¬ 
tain harmony of the terminology with the planar slab 
example. In the nonactive scenario, as before for the pla¬ 
nar case, the normal (radial) permittivity and permeabil¬ 


ity are assumed to be real 


(s' , 

\^rp' 


= Re 


a-jb 


a2+62 


while the other material parameters are given by (23). 


In the PML route, we again assume that the normal 
components of permittivity and permeability are slightly 
more lossy than it is dictated by the ideal PML rule 
(given in ( [2^ ) in order to avoid possible numerical in¬ 
stabilities owed to the limiting nature of the ideal CML 

(4p' = 


B. Performance of the Proposed Structure 


In Figs.[^and[^ we consider the case of two cylinders 
having material parameters corresponding to the three 
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(c) 


FIG. 12. The spatial variation of the axial magnetic field 
Hz{x,y) on the map {x,y) for the three considered DPS 
routes: (a) the isotropic route {e'^p — a — jb), (b) the nonac¬ 
tive route [s'rp = Re [l/(a — jb)] — a/{a^ + 5^)), and (c) the 
PML route {e'rp — l/(a — jb) — jS). Plot parameters: a = 2, 
kor = 20, kod = 3, 6> = 90°, 5 = 0.05, S = 0.005, and n = 1. 



(c) 


FIG. 13. The spatial variation of the axial magnetic field 
Hz{x,y) on the map {x,y) for the three considered DNG 
routes: (a) the isotropic route = a—jb), (b) the nonactive 
route {e'rp — Re [l/(a — jb)] — a/{a^ + ^‘^))^ (c) the PML 

route (^e'rp = l/{a — jb) — jS) . Plot parameters: a = —2, 
kor = 20, kod = 3, 6> = 90°, b = 0.05, S = 0.005, and n = 1. 


routes described in the previous subsection (“isot ropic”, 
nonactive, and PML routes) and depicted in Fig. |ll(b)| 
We position the source line symmetrically in-between two 
cylinders of radius r which are separated by distance d. 
When the source is that depicted in Fig. namely a 
TM dipole line source, we again choose 0 = 90° for m ax- 
imal excitation of the system. In Figs. 12(a) - p!2(c) 


we 


consi der the DPS (PML) cases with a = 2, while in 
Figs. |13(a)f T3(^ the DNG (CML) ones with a = —2. 
By inspection of the contours, one directly infers that 
the DNG approach is much more efficient in terms of 
absorption: in the DNG case, practically all the power 
created by the source is absorbed by the cylinders. On 
the contrary, in all the DPS cases a considerable amount 
of power propagates away from the cylinders. Similar 
effect was observed for the DNG cylinders even if the 
dipole line excitation was replaced by a 2 :-directed cur¬ 
rent line, despite the fact that the incident field is nonzero 
along the y axis. The difference between Figs. |12(b)| 
and between Figs. |13(b)|13(c)| is negligible due to 


small losses b com bined with the choice of fairly large \a\. 
Figs. |13(b)|13(c)| show that surface plasmon modes are 
strongly excited due to smaller overall losses compared 
to the other cases. In the nonactive and PML routes 
with a = —2 shown in Figs. |13(b)| and |13(c)[ surface 
plasmon modes are strongly excited on the surface of the 
cylinders which d oes not happen for the isotropic DNG 
case shown in Fig. |13(a)] dne to higher losses in the radial 
material parameters compared to the other cases. 

In Fig. 14(a)| we show the variations of the absorbed 


power P normalized by Pq, which corresponds to the ab¬ 
sorbed power in the DPS-PML case with a = 1, with 
respect to a. Given the fact that the conventional PML 
absorbs only propagating waves, the physical meaning of 
Po is identical to the corresponding one of the planar slab 
cases. Regardless of the considered route (isotropic, non¬ 
active, PML), there is a substantial switch in absorption 
efficiency from a < 0 (DNG case) to a > 0 (DPS case): 
the absorbed power P is always much higher for a < 0. 
For positive values of a, all routes lead to similar results. 
The higher absorption in the isotropic shell case is simply 
due to the lossy radial component (as opposed to non¬ 
active/active Srp used in the other two cases). When it 
comes to the isotropic shell case, the maximal absorp¬ 
tion is achieved when a = —1, since the power drops 
quickly as a deviates from —1 in either direction; similar 
behavior has been observed also for the planar slab case, 
as well as for the double-negative spher^^. It should 
be noted that at a = ±1 practically the only difference 
between the material parameters given by the three sce¬ 
narios are the losses in the radial permittivity and per¬ 
meability. Furthermore, it should be stressed that in the 
nonactive DNG case, the maximum absorbed power (ex¬ 
hibited at a = —2.5) is smaller than in the PML scenario 
but larger than in the isotropic shell case, as could be 
expected. Similarly to the planar case , the h ighest P is 
recorded for the PML route. In Fig. |14(b)| we present 
the escaped power P+ with respect to the real part of 
the transversal relative permittivity a. P+ denotes the 
power carried away by the surviving fields at p' +oo. 
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(b) 

FIG. 14. (a) The normalized absorbed power P/Pq by the two 
cylinders and (b) the escaping power P^/Pq from the system 
as functions of the material parameter a in the three differ¬ 
ent scenarios (isotropic, nonactive, PML). Plot parameters: 
kor = 20, kod = 3,0 = 90°, b = 0.05, 6 = 0.005. 


The values of P+ for a < 0 are greatly diminished com¬ 
pared to the DPS cases, a feature that demonstrates the 
higher absorption effectiveness of the DNG structures. 
For \a\ > 1, the escaped power is the smallest for the 
isotropic scenario which is, again, due to higher losses 
compared to the other cases. 

Finally, in Figs. we examine absorption effective¬ 
ness of double-negative conjugate matched cylinders il¬ 
luminated by a single propagating plane wave. Accord¬ 
ing to the results of Ref. [191 fhe limiting case of the 
ideal CML, the effective absorption width can be arbi¬ 
trarily large, although the sources of illumination are in¬ 
finitely far from the cylinder. We consider only the PML 
route for a = 2 (DPS, Fig. 15(a)[ ) and a = —2 (DNG, 
Fig. 115(1^ . Note that in both cases, there are practically 
no reflections from the cylinder and a huge shadow is 
formed behind. Furthermore, resonant plasmonic modes 
are excited on the surface of the cylinder in the DNG 
case, which increases the shadow due to stronger field 
coupling. Although the excitation is a single propagat¬ 



FIG. 15. Plane-wave illumination. The spatial variation of 
the axial magnetic field Hz{x,y) on the map {x,y) for the 
PML route [e'rp = l/{a — jh) — jS) with (a) a = 2 and (b) 
a = —2. Plot parameters: kor = 20, kod = 3, b = 0.05, 
S = 0.005, and n = 1. 


ing plane wave, coupling to predominantly evanescent 
resonant modes of the absorbing cylinder is possible be¬ 
cause the surface is curved. In this particular example, 
the power absorbed by the cylinder is about 18% higher 
in the DNG case compared to the DPS case (the latter 
case corresponds to the ideal black body as defined in 
Ref.^. In other words, even in the absence of incident 
evanescent plane-wave modes we can achieve improved 
absorption using DNG cylinders as compared to the per¬ 
fectly non-reflecting DPS cylinder, thanks to the afore¬ 
mentioned activation of resonant surface plasmon modes. 
However, for the same material parameters, the effect is 
quite modest compared to the increased power absorp¬ 
tion when evanescent wave modes are strongly present in 
the spatial spectrum of the source. 


V. PRACTICAL REALIZATION OF THE CML 
CONCEPT 

As is clear from the above results, two scenarios are 
promising for realizations of the proposed conjugate 
matched bodies: the nonactive route and the PML route. 
In the nonactive route, the required material parame- 
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ters correspond to well-studied uniaxial double-negative 
(backward-wave) uietauiaterial^^. A wide variety of 
micro- and nanostructures which realize double-negative 
response in different frequency ranges have been pro¬ 
posed and studied in the literature. For the new ap¬ 
plication proposed here, the materials should have a con¬ 
trollable degree of anisotropy: in particular, the case 
when the absolute values of the tangential components 
are much larger than those of the normal components is 
of interest. With this in view, one of the most promising 
topologies is the so-call fishnet structurd^^ll^. This is a 
multilayer structure where thin metal sheets with peri¬ 
odically positioned holes (usually of the square or round 
shape) are separated by dielectric spacers. 

While the response in the tangential plane is strongly 
modified by large induced currents flowing along the 
metal sheets, the response to fields along the normal di¬ 
rection is modified only due to quasi-static electric po¬ 
larization of thin sheets. In the microwave domain, the 
normal components of the material parameters usually 
have positive real parts. To achieve the necessary neg¬ 
ative permittivity value of the normal component, it is 
possible to insert an array of thin metal wires passing 
through the holes in the fishnet layers. If the wires are 
made of a good conductor, the microwave response will 
correspond to a negative and nearly real value of the 
normal component of the permittivity tensor, exactly as 
required for the CML realization. 

At infrared and optical frequencies, it is possible to 
exploit the negative permittivity of metals from which 
the fishnet layers can be fabricated. In this case, metal 
wires are not needed, and the values of the normal com¬ 
ponent of the permittivity can be controlled by choos¬ 
ing the thickness ratio of the metal and dielectric lay¬ 
ers in the fishnet structure. In particular, if rectangu¬ 
lar holes are etched in a layer of Magnesium fluoride 
(MgF 2 ) sandwiched between silver slabs, one can achieve 
negative refractive inde:xP^ at the optical frequencies un¬ 
der the necessary homogenization conditions^. Further¬ 
more, nonlinearity of liquid crystals can be exploited in 
order to develop tunable fishnets metamaterials whose 
reflection/transmission (and accordingly their effective 
^rtii^rt) is externally controlled. More generally speak¬ 
ing, it is well known that if we insert metallic particles 
such as cylindrical pins, strips, spirals or flakes, which are 
having negative permittivity in the visible spectrunP^, 
we caiP^ obtain materials with Re[5rn] < 0- 

The conclusion that the performance of nonactive- 
route CML improves when the absolute value of the real 
part of the tangential permittivity increases (Section |lll| ) 
is very important for choosing practically realizable ma¬ 
terial parameters. As was discussed above, high absorp¬ 
tion of energy from evanescent fields requires that the 
overall losses in the medium are small (this is in fact 
clear already from eq. (§)• Realization of passive low- 
loss DNG metamaterials is a chall eng e, although some 
successful approaches are known, e.g.®. However, for the 
nonactive PML with large magnitudes of negative a, high 


absorption is observed even when the imaginary part of 
the tangential permittivity b is not very small and corre¬ 
sponds to practically realizable structures. The normal 
component of the permittivity in this case corresponds 
to epsilon-near-zero m ateria ls, which are also known to 
be practically realizabl d^^ ^ ^^ l 

The required overall small level of losses can achieved 
also using active (pumped) structures. This can be 
done by using active parts such as optically pumped 
laser dyes to balance the losses of the background 
structurd^. In this way, one can realize an effective DNG 
medium possessing additional negative normal compo¬ 
nents (Re[£^^] < 0) with controllable loss factors of the 
tangential components: Im[5rt],Im[/i^t] < 0- Active con¬ 
trol of the normal permittivity and permeability com¬ 
ponents can be used to implement also the active PML 
route towards GML. Based on the fishnet-wire medium 
approach described above, instead of metallic rods or spi¬ 
rals one can insert nano-generatord^ or pumpd^into the 
holes of the background structure which can provide en¬ 
ergy to the system and make the normal component of 
the effective parameters active. 


VI. CONCLUSIONS 

In this paper, we have introduced the concept of 
double-negative conjugate matched layer (CML), which 
has the property of acting as the ideal sink for energy of 
electromagnetic fields. Arbitrary propagating waves illu¬ 
minating such layer produce no reflections and are fully 
absorbed in the medium. Evanescent fields, existing in 
the vicinity of small sources or scatterers, induce very 
strong resonant fields, oscillating in the vicinity of the 
CML surface. The energy stored in the reactive evanes¬ 
cent modes sinks into the layer and is transformed into 
heat or delivered to loads in its constitutive elements. 

There is an analogy of this phenomenon and the 
circuit-theory concept of the ideal voltage or current 
source. An ideal voltage source delivers infinite power 
to the load in the limit of infinitely small resistive load 
[eq. (©], because the current through a resistor tends to 
infinity when the resistance tends to zero. Likewise, for 
every evanescent mode the sum of the wave impedances 
of fields in free space and that in the CML is real and 
can be made very small. Thus, the layer dramatically en¬ 
hances the near fields created by given sources and sacks 
their energy into the material layer. This phenomenon 
can be compared also with the Purcell effect of enhance¬ 
ment of spontaneous emission rate of small sources due to 
the presence of resonant a nd eff iciently radiating bodies 
in the vicinity of the sourcd^^EZl. However, resonant cav¬ 
ities or antennas used to enhance radiation from a small 
emitter work only at one or a few resonant modes of the 
resonator. The introduced CML concept brings the Pur¬ 
cell effect to the limit of extracting all available energy 
from a small source. In this paradigm, every mode of the 
source field resonates with a corresponding mode of the 
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CML, so that the complete system of the source and the 
absorbing body is tuned for the optimal power extraction 
from all modes of the fields created by the source. 

Although the ideal double-negative conjugate matched 
layer re quir es active inclusions for its practical 
realizatiorP^ (the normal components of the material 
parameters have the imaginary parts corresponding to 
active media), the flexibility of the CML concept sug¬ 
gests a more practical possibility for realizations, where 
only passive materials are used: the nonactive route, 
described in this paper, requires double-negative lossy 
values of the tangential components and low-loss neg¬ 
ative values of the normal components of the permit¬ 
tivity and permeability matrices. The results of exten¬ 
sive studies of various possibilities to design and realize 
double-negative (backward-wave) materials can be used 
to realize the proposed structures. While the earlier in¬ 
troduced isotropic double-negative desigrP^ requires very 
small levels of losses, which are probably not possible to 
achieve in practice, at least with the use of only passive 
materials, the radially uniaxial double-negative structure 
proposed here appears to offer possibilities to overcome 
this limitation. We hope that the results of this study 
will have interesting and important implications for a 
wide variety of applications: in light harvesting, thermal 
emission control, heating and cooling devices, stealth, de¬ 
coupling of radiators, receiving and transmitting anten¬ 
nas, and other microwave and optical devices. 


written for this 2D problem as 


Ao,inc(^a:5 ^y) 


jkoijopie 

Z-2 I iL2 _ 7.2 ’ 

I '^y ^0 


(A.l) 


where is the electric dipole moment linear 

density. Taking the inverse Fourier transform 

(A).inc(a;, ky) = ^ Ao,i„c(fca:, ky)e~^'^^^dka^ with 

respect to we obtain 




^ Jo 


(A.2) 


where /3o = Po{^y) = —j When ky = 0, (A.2) 

reduces to the well-known expression for the vector po¬ 
tential of a uniform electric dipole moment sheet. The 
magnetic field of the same line of dipoles is found using: 


"^0,inc(^5 ^y) — V X wAo^inc (^: ) 5 

Mo 


(A.3) 
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Appendix: Expression for the Incident Field Created 
by the Dipole Line 

In the spectral domain, the free-space vector potential 
created by the line of dipoles shown in Fig. can be 


where V = :k-^—yjky (recall that there is no dependence 
on z). The result reads 


"^0,inc(^: ) 


2 


~ k 

cos 0 — sin 0 
_Po 


(A.4) 

where pi = |p/|. Because the vector pi lies in the xy 
plane, the incident magnetic field has just a single com¬ 
ponent along the 2 : axis, i.e., the incident waves are TM- 
polarized plane waves. 
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